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ABSTRACT 


We discuss how laboratory experiments can be used to place constraints 
on possible variations of the fine structure constant a in the observationally 
relevant redshift interval z ~ 0 — 5, within a rather general theoretical 
framework. We find a worst case upper limit for A a/a of 8 x ICC 6 for z < 5 
and A a/a of 0.9 x 10 -6 for z < 1.6. The derived limits are at variance with 
the recent findings by Webb et al. (1998), who claim an observed variation of 
A a/a = —2.6 ± 0.4 x 10~ 5 at 1 < z < 1.6. 

subject headings : cosmology: theory - physical data and theory: relativity, nuclear 
reactions, nucleosynthesis, abundances 
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1. Introduction 

In an interesting paper, Webb et al. (1998) describe an observational method for 
investigating possible time and space variations in the fine structure constant a, using 
quasar absorption systems (see also Drinkwater et al. 1998). Webb et al. present intriguing 
evidence suggesting that a was smaller in the past (for redshifts z > 1). The claimed 
fractional change is (for z > 1), A a/a = —2.6 ± 0.4 x 10“ 5 . 

Since a varying value of a can, in principle at least, have significant consequences 
for cosmology (e.g. concerning recombination), we examine in the present letter possible 
variations in a in a broader context, although still under a particular set of assumptions. 
We also show that in the context of the theory presented in this paper the result obtained 
by Webb et al. appears to be in conflict with results of tests of the equivalence principle. 


2. By How Much Can a Change? 

The experimental constraints on variations of a were explored extensively by Bekenstein 
(1982) by constructing a dynamical theory relating variations of a to the electromagnetic 
fraction of the mass density in the universe. Bekenstein’s (1982) framework for a variability 
was based on very general assumptions: covariance, gauge invariance, causality, time 
reversal of electromagnetism, and that gravitation is described by a metric of spacetime 
which satisfies Einstein’s equations. He obtained the following equation for the temporal 
variation of a (adopting a Robert son-Walker metric; we correct here a misprint in the 
original paper) 

(a 3 e/ej = -a 3 ((l 2 /hc)p m c 4 . (1) 

Here e = (a/atoday) 1 ^ 2 , l is a length scale of the theory, p m is the total rest mass density of 
matter, a is the expansion scale factor, and ( = (m ntem )/m p is a dimensionless parameter 
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which measures the fraction of mass in the form of Coulomb energy of an average nucleon, 
(■ m n ,em ), compared to the free proton mass, m p . In order to be able to integrate eq. (1), 
Bekenstein (1982) assumed that ( (which is affected by transformations of hydrogen into 
heavier nuclei in stars) is constant. He was able to show that this assumption is reasonable 
up to redshifts z 7S 1, by the following consideration. The proton conversion rate (via the 
reaction 4H 1 —» He 4 ) can be estimated from the mean luminosity density (e.g. Davis, Geller 
& Huchra 1978) to be ~ 8 x ICC 28 cm -3 s -1 . Since a minimal nucleon number density in the 
universe can be determined from nucleosynthesis to be (e.g. Olive, Schramm & Steigman 
1981) ~ 10 -7 cm -3 , the timescale for proton conversion (and thereby for a change in Q) is 
T ;> 10 20 S 

' conv ~ 0 • 

The timescale for changes in p m or a is rsj H- 1 10 17 s. Thus, for as far back as the 

luminosity of galaxies does not change significantly (up to z 7$ 1), the assumption of a 
constant ( is a reasonable one. 

However, it is interesting to push this analysis to higher redshifts (e.g. the main change 
in a claimed by Webb et al. (1998) occurs at z ^ 1) by extending Bekenstein’s (1982) 
analysis with a determination of the behavior of ( at higher redshifts. We achieve this as 
follows. As explained above, the rate of change of the Coulomb contribution to the baryon 
mass depends on the rate of conversion of hydrogen into helium. Following Bekenstein we 
find: 

C ~ 1.2 x 10" 2 (X + -Y) , (2) 

3 

where X and Y are the mass fractions in hydrogen and helium, respectively. In order to 
determine the conversion rate we adopted three possible star formation rate (SFR) histories 
for the universe: (i) one that follows the Madau et al. (1996) curve, (ii) one that rises from 
the present up to a peak at z ~ 1 and then is flat at the peak value up to z ~ 5 (as is 
perhaps suggested by the COBE Diffuse Infrared Background experiment; Hauser et al. 
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1998; Calzetti & Heckman 1998), and (iii) a SFR that is constant (at the peak value) for 
all redshifts (as a limiting case). These star formation histories were used to synthesize 
a mean stellar population of the universe as a function of time, using the Bruzual and 
Chariot (1993) evolutionary models. The bolometric luminosity from the mean stellar 
population was then used to calculate the hydrogen to helium conversion rate. In Fig. 1, 
we show the proton conversion rate as a function of redshift up to z = 5, for the three 
assumed SFR histories. As can be seen from the figure, the timescale for proton conversion 
satisfies T conv ^ 10 19 s up to z = 5. Since r conv H ~ 1 , we can still take ( to be constant in 
integrating eq. (1). Indeed direct calculation of eq. (H) shows that it varies by less than 1 % 
over the redshift interval that we consider. Thus we obtain: 


e/e = -C (l' 2 c 3 /h)p m (t - t c ) 


( 3 ) 


where t c is an (unknown) integration constant. From eq. @ with A" = 0.74 and Y = 0.26 
we find for ( the value 1.3 x 1CU 2 (see also Bekenstein 1982). We will consider in the 
following the case of a low (fi ~ 0.2, see below) Universe. The results however would 
not change significantly in a = 1 Universe. Eq. (^j) can be integrated in time by taking 
into account that for z << 4(f2 _1 — 1) one has p m ~ Q m0 p c (t 0 /t) 3 , where Q mo is the baryon 
fraction at the present time t 0 , and p c is the critical density p c = 3Hq/(8ttG), and that 
\e — 1| << 1. We find: 


e(z) — 1 ~ 


,3C 

87T 


l 


( 1 + 


t r t\ 


c L 0 


to 


2 t 2 (z) t(z) 2 t 0/ 


( 4 ) 


where L p = (Gfr/c 3 ) 1 / 2 is the Planck length. We have verified that the above analytical 
approximation is applicable to the redshift interval of interest. As an example, in Fig. 2 we 
plot for the case t c = to the difference between the expression of eq. (4) and the result of a 
direct numerical integration of eq. (3). Eq. (]4|) has two quite different regimes depending 
on the value of t c /t Q . Let us consider first the case \t c \ <t Q which is perhaps the most 
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physically founded since there is no a priori reason for having \t c \ » t Q . We find: 

|e(z < 5) - 1| < 6.2 x 10” 3 f2 mo (-pj . (5) 

In order to obtain an absolute upper limit on the variability, we will assume that all the 
matter density is nucleonic, by taking = Vt mo ~ 0.2 (e.g. Garnavich et al. 1998; Perlmutter 
et al. 1998). We can also use the results of the Eotvos-Dicke-Braginsky experiments 
(Eotvos, Parker & Fekete 1922; Roll, Krotkov & Dicke 1964; Braginsky & Parnov 1972), 
designed to test the equivalence principle, to infer 1/L p £ 10” 3 . Thus, our final upper limit 
is: 

|e(-< 5) - 1| < 1.2 x 10“ 9 . (6) 

If we assume instead that \t c \ >> t 0 we can set an upper limit to \t c \ by considering 
constraints on \e/e\. From Shlyakhter’s (1976) analysis of the Oklo natural reactor we have 
|e/e| < 0.5 x 10 ~ 7 H a . More recently, Damour and Dyson (1996) have derived the less strong 
but more robust limit of |e/e| < 3.4 x 10~‘H 0 , which we will adopt in the following. When 
used in combination with eq. (j3j) and with Q mo = 0.2 this gives: 

\(l/L p ) 2 t c H 0 \ < 6.6 x 10~ 3 . (7) 

The term involving t c in Eq. ([!]) becomes dominant and we obtain: 

\e( z < 5) — 1| < 8 x 10“ 6 . (8) 

The above limits become stronger as the redshift decreases. For example, at the 
redshift of interest for the Webb et al. (1998) result (z=1.6) one would find instead 
| e(z = 1.6) — 1| < 1.6 x 10 -1 ° if \t c \/t a < 1 or \e( z = 1-6) — 11 < 9 x 10”' otherwise. Both 
limits, regardless of the value of t c , are at variance with the Webb et al. (1998) result. 
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3. Discussion 

The analysis of experimental constraints for the variation of the fine structure constant 
within the framework of rather general dynamical theories provides upper limits to any 
change in a which are at variance with the Webb et. al (1998) detection. We have verified 
that varying the cosmological parameters (i.e. the values of fi, Hq, primordial helium 
fraction) does not affect our constraints significantly. 

We should note that variations in a in the context of the theory discussed in the 
present paper result in no modifications to the Planckian spectrum of black body radiation. 
In fact one could think of varying a as a varying permittivity of the vacuum. This would 
simply be equivalent to changing the expansion factor a(t) to a different function, which 
cannot modify the Planck spectrum. 

Variations of the fine structure constant a could be driven also by a different dynamics. 
In theories with extra dimensions (e.g., Kaluza-Klein or superstrings theories) variations 
in the observed value of a can be induced by variations in the length scale of the compact 
dimensions (see, e.g., Marciano 1984, Barrow 1987). The effects of extra dimensions were 
discussed in detail by Barrow (1987) who found upper limits to \a/a\ of < 10~ 8 H 0 for 
both Kaluza-Klein and superstring theories. This upper limit is also much smaller than the 
claimed detection by Webb et. al (1998). However, note that in these theories the variations 
of all physical constants are inter-related. 

Finally, we would like to mention that the calculations leading to eq. (f|) would have 
to be modified if a significant component of the dark matter in the universe were to decay 
electromagnetically at redshift z ^ 1 (e.g. Sciama 1982, 1997). We hope that the present 
letter will inspire more attempts to determine possible variations of a observationally. 

ML acknowledges support from NASA Grant NAG5-6857. We thank Yacov Bekenstein 


for very useful discussions and an anonymous referee for a careful reading of the manuscript. 
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Fig. 1.— We plot the proton conversion rate as a function of redshift for 
model (solid line), the Madau et. al cosmic SFR history (dotted line) and 
rising up to z — 1 and then staying constant at the peak value (dashed line). 


constant SF 
SFR history 


Fig. 2.— Difference between the analytical approximation and the numerically integrated 
value of e. For redshift z < 10 the difference between the two is negligible. 







